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■ Abstract 

i^h ' We provide a priori error estimates for the spectral and pseudospectral 

C$ , Fourier (also called planewave) discretizations of the periodic Thomas-Fcrmi- 

von Weizsacker (TFW) model and of the Kohn-Sham model, within the local 
density approximation (LDA). These models allow to compute approximations 
of the ground state energy and density of molecular systems in the condensed 
phase. The TFW model is stricly convex with respect to the electronic density, 
^ | and allows for a comprehensive analysis (Part I). This is not the case for 

. the Kohn-Sham LDA model, for which the uniqueness of the ground state 

electronic density is not guaranteed. Under a coercivity assumption on the 
second order optimality condition, we prove in Part II that for large enough 
energy cut-offs, the discretized Kohn-Sham LDA problem has a minimizer in 
0^ ' the vicinity of any Kohn-Sham ground state, and that this minimizer is unique 

, up to unitary transform. We then derive optimal a priori error estimates for 

■ both the spectral and the pseudospectral discretization methods. 



1 Introduction 

Density Functional Theory (DFT) is a powerful method for computing ground state 
electronic energies and densities in quantum chemistry, materials science, molecular 
biology and nanosciences. The models originating from DFT can be classified into 
two categories: the orbital-free models and the Kohn-Sham models. The Thomas- 
Fermi- von Weizsacker (TFW) model falls into the first category. It is not very much 
used in practice, but is interesting from a mathematical viewpoint. It indeed serves 
as a toy model for the analysis of the more complex electronic structure models 
routinely used by Physicists and Chemists. At the other extremity of the spectrum, 
the Kohn-Sham models are among the most widely used models in Physics and 
Chemistry, but are much more difficult to deal with. We focus here on the numer- 
ical analysis of the TFW model on the one hand, and of the Kohn-Sham model, 
within the local density approximation (LDA), on the other hand. More precisely, 
we are interested in the pseudospectral Fourier, more commonly called planewave, 
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discretization of the periodic version of these two models. In this context, the simu- 
lation domain, sometimes rcfercd to as the superccll, is the unit cell of some periodic 
lattice of R 3 . In the TFW model, periodic boundary conditions (PBC) are imposed 
to the density; in the Kohn-Sham framework, they arc imposed to the Kohn-Sham 
orbitals (Born-von Karman PBC). Imposing PBC at the boundary of the simula- 
tion cell is a standard method to compute condensed phase properties with a limited 
number of atoms in the simulation cell, hence at a moderate computational cost. 

This article is organized as follows. In Section 2, we briefly introduce the functional 
setting used in the formulation and the analysis of the planewavc discretization 
of orbital-free and Kohn-Sham models. In Section 3, we provide a priori error 
estimates for the planewave discretization of the TFW model. Our estimates refine 
and complement some of the results given in [10]. In Part II, we deal with the 
Kohn-Sham LDA model. 



2 Basic Fourier analysis for planewave discretiza- 
tion methods 

Throughout this article, we denote by V the simulation cell, by TZ the periodic lattice, 
and by TZ* the dual lattice. For simplicity, we assume that T = [0, L) 3 (L > 0), but 
our arguments can be straightforwardly extended to rectangular simulation cells 
(r = [0,L X ) x [0,L y ) x [0,L 2 )). For T = [0,L) 3 , TZ is the cubic lattice LI 3 , and 
1Z* = ?fl 3 . For k G TZ*, we denote by e k (x) = \T\- 1 / 2 e ik ' x the planewave with 
wavevector k. The family (ek)keK* forms an orthonormal basis of 

L%(T,C) ■= {u 6 Lf oc {m 3 ,C) | u ^-periodic} , 

and for all u G £|(r,C), 



(x)= ^u k e k {x) with u k = (e k , u) L 2 # = \T\ 1/2 / u(x)e ik ' x 
ken* Jv 



dx. 



In our analysis, we will only consider real valued functions. We therefore introduce 
the Sobolev spaces of real valued functions 

H S # {T) := J u{x) = Y,u k e k (x)\ ^ (1 + |fc| 2 ) s |w fe | 2 < co and Vfc, c_ fc = c% \ , 
I ken* ken* ) 

s£l, endowed with the inner products 

ken* 

For N c G N, we denote by 



V N ,. = { y~] c k e k | Vfc, c_ fc = c* k ) (1) 




(the constraints C- k = c| imply that the functions of Vn c are real valued). For all 
s£l, and each v G H^(T), the best approximation of v in Vn c for any TJ^-norm, 
r < s, is 

ITv c v = 2J " fcefc - 

ken* | \k\<zg-N c 
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The more regular v (the regularity being measured in terms of the Sobolev norms 
H r ), the faster the convergence of this truncated series to v: for all real numbers r 
and s with r < s, we have for each v £ H^(T), 



\v-U Na v\\ Hi = min \\v-v Ne \\ H r < — N~^\\v -U Nc v\\ H ^ 



< [^) N-^\\v\\ Hk (2) 

For N g £ N \ {0}, we denote by FFT ' Jv s the discrete Fourier transform on the 
carterisan grid Qn 9 '■= jf- Z 3 of the function <t> £ C^iT). Recall that if <j) = 

12ken* < t ) g e g e the discrete Fourier transform of <f) is the iV s 7£*-periodic 

2FFT N f7PFT,N g ^ , 

sequence (p rr± '"n = (0 fe )ken* where 

r" 3 = ^ E - i r r 1/2 E W- 

9 xtg Ng nr Ken* 
We now introduce the subspaces 



W£ = 



Span Uv \l£^Z, \l\ < ^ f^^) } (N g odd), 

Span L ily | J G ^Z, \l\ < ^ | 8 C(e i%N ° y/L + e -^ N ^ L ) (N g even), 

(W}P g £ C£([0, L)) and dim(W^) = N g ), and = ® ® Note 
that Wjy D is a subspace of H^(T) of dimension TV 3 , for all s £ R, and that if iV g is 
odd, 

= Span |e fc | fc £ TV = ^Z 3 , (fe^ < ^ (^^) } odd )- 

It is then possible to define the interpolation projector Xjv s from C^(r) onto WfP 
by [In ((j))](x) = </>(x) for all x £ Gn s - It holds 

V0ec°(r), fi Ng (<f>) = E (at)^)- ( 3 ) 

The coefficients of the expansion of Xjv (4>) in the canonical basis of WfP is given 
by the discrete Fourier transform of <f>. In particular, when N g is odd, we have the 
simple relation 

^)Hlf /2 E ?™ 9 e fc (A^odd). 

ken* \\k\^<^{^-) 

It is easy to check that if (j> is real- valued, then so is Xjv 9 (</>). 

We will assume in the sequel that iV fl > AN C + 1. We will then have for all i>4at c G 



/ U 4W C = E \Tr) V ^A x ) = [ X N, 

Jr ,ee«„nr V^/ Jr 



(V4N C ). (4) 



The following lemma gathers some technical results which will be useful for the 
numerical analysis of the planewave discretization of orbital-free and Kohn-Sham 
models. 
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Lemma 2.1 Let N c £ N* and N a £ N* such that N a > AN C + 1. 



i. Le£ be a real-valued function of C^.(T) and Vn c and wn c be two functions 
of Vn c ■ Then 



lN g (Vv Nc w Nc ) = / l Ng (V)v Nc w Nc 
r Jr 



lN g {V\v Nc \ 2 ) 



< \\V\\l~\\v Nc \\] 



(5) 
(6) 



2. Let s > 3/2, < r < s, and V a function of H^(T). Then, 
11(1-2^)^)11^ < C r , s N r g ~ s \\V\\ Hi 



\u 2Nc (i N (v))\\ L2 < / i Ng (\v\ 2 ) 



1/2 



\TL 2Ne (l Ng (V))\\ He < (l + C s , a )\\V\\ H s 



(7) 

(8) 
(9) 



for constants C r>s independent ofV. Besides if there exists m > 3 and C £ R+ 
such that \Vu\ < C|fc| _m , then there exists a constant Cy independent of N c 



and N g such that 



\lL 2Nc (l-I Ng ){V)\\ Hr < CvK+^N^ 1 



(10) 



3. Let <j) is be a Borel function from R + to R such that there exists £ R+ for 
which \<j>(t)\ < C (l + t 2 ) for all t £ R+. Then, for all v Nc £ V Nc , 



TnMIvnJ 2 )) 



< C A [\T 



\VN C \\% 



(11) 



Proof For z 2 n c £ V2at c , it holds 
I Ng (Vz 2Nc ) 



E (lN g (V))(x)z 2Nc (a 



T Ng (V) z 2Nc 



(12) 



smce I N g (V)z 2 N c £ VW C - The function vn c wn c being in V 2 n c , (5) is proved. More- 
over, as IvnJ 2 £ ViN c , it follows from (4) that 



lN g {V\v Nc 



E uH 



< ||V|| 



L x 



E (w ) '^ (a;)|2 
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Hence (6). The estimate (7) is proved in [4]. To prove (8), we notice that 
HIWI^F))!!^ < \\lN g (V)f L2 

w 

[iN a {v)y{i Ng {v)) 

E (lN a (V)){x)*{I Ng {V)){x) 

E i^)i 2 

xeg Ng nr 
In 9 (\V\ 2 ). 

The bound (9) is a straightforward consequence of (7): 

\\H 2Nc (lN 9 (V))\\ Hi < \\lN g (V)\\H^<\\V\^ + \\(l-I Ng )(V)\\H^<(l + C s ^\\V\\ H ^ 

Now, we notice that 

k<=1f | \k\<ig-N a 



E ( E wl 

ken* I \k\<^N r \Keii' 



(13) 



From (13), we obtain 



\U 2Nc (1-Xn 3 )(V)\\ hb 



E ( i + \ k \ 2 r 

fees* | |fe|<^Ar c 



E 

KEK*\{0} 



V (l + |fc| 2 )^| max 

f—f . / fc€R* I \k\<i^N € 

^k£W | \k\<i£N a 



E % 

K£K*\{0} 



+N g K 



On the one hand. 



E (i 



|2\s 



32tt /47T 



jv c -»oo 2s + 3 V L 



2 s 



N 2s+3 



and on the other hand, we have for each k 6 TZ* such that |fc| < ^N c , 



E Vk+NgK 

KeK'\{0} 



z c E 

ife-R.*\{o} 



|fc + N g KY 



where 



C 



max N 

y \<i/2 ^ \y-K\ r 



The estimate (10) then easily follows. Let us finally prove (11). Using (3) and (4), 
we have 



1nM\vn\ 2 )) 



xe-^z 3 nr 



E (w) ^M*)! 2 ) 



L 



c y r (i + |^| 4 ) = c (|r| + ||^||^) 



This completes the proof of Lemma 2.1. 



□ 



3 Thomas-Fermi-von-Weizsacker model 



In the TFW model, as well as in any orbital-free model, the ground state electronic 
density of the system is obtained by minimizing an explicit functional of the density. 
Denoting by Af the number of electrons in the simulation cell and by 

3W = |p>o|v?e##(r), f p = j A 

the set of admissible densities, the TFW problem reads 

7 TFW = mf{£ TFW (p), pe^}, (14) 

where 



Cw is a positive real number (Cw = 1, 1/5 or 1/9 depending on the context [5]), 
and Ctf is the Thomas- Fermi constant: Ctf = -y(37r 2 ) 2 / 3 . The last term of the 
TFW energy models the periodic Coulomb energy: for p and p' in i?^ 1 (r), 

feG7J*\{0} 

We finally make the assumption that V lon is a periodic potential such that 

3m > 3, C > s.t. Vfc G K*, \V^ on \ < C\k\- m . (15) 

Note that this implies that V" lon is in H m - 3 / 2 - e (T) for all e > 0. It is convenient to 
reformulate the TFW model in terms of v = ^fp. It can be seen that 

I TFW = inf |b tfw («), v G ff* (T), jf | V | 2 = AA| (16) 

where 

£ TFW H = y/ r N 2 + c TF | r H 10 / 3 + | r ^H 2 + ^ r (H 2 ,H 2 ). 



6 



It is well known [6] that (14) has a unique minimizer p°, and that the minimizers 
of (16) are u and —u, where u = ■ Besides, the function u is in 
for any e > (and therefore in Cjt(r) since to + 1/2 — e > 7/2 for e small enough), 
is positive everywhere in T and satisfies the Eulcr equation 

-^Au + ( lc TF u 4/3 + V ion + V3° ul0mh ) u = Xu 



2 — 
for some Ael, where 

V^° ulomb (z) = 4?r ]T \k\- 2 p k e k (x) 
keiz*\{o} 

is the periodic Coulomb potential generated by the periodic charge distribution p. 
Recall that v^ Coulomb can also be defined as the unique solution in H^(T) to 



. A y -Coulomb = 4?r ^ _ |p|-l J 



^Coulomb = Q 

r 

The planewave discretization of the TFW model is obtained by choosing 

1. an energy cut-off E c > or, equivalently, a finite dimensional Fourier space 
Vn c , the integer N c being related to E c through the relation N c := [y / 2E c L/2n]: 

2. a cartesian grid Qn with step size L/N g where N g g N* is such that N g > 
4iV c + 1, 

and by considering the finite dimensional minimization problem 

a 



I™ g = inf {e™(v Nc ), V Nc G y Wc , jf |^J 2 = AT} , (17) 



where 



^f W KJ = ^ / |V^J 2 + C TF / J^d^l^) + / JjVg(F io n)| ^| 2 

r *f r */ r 

+iiMKJ 2 > J vJ 2 ), 

Xjv s denoting the interpolation operator introduced in the previous section. The 
Euler equation associated with (17) can be written as a nonlinear eigenvalue problem 

Vv Nc € V Nc , ({H" g Nc N ^ 2 U NctNg - XN c ,N g )UN c ,N g ,V Nc } H -^ H i # = 

where we have denoted by 



lomb 



the pseudo-spectral TFW Hamiltonian associated with the density p, and by Xn c ,n 
the Lagrange multiplier of the constraint J r |vjvj 2 = M. We therefore have 



~-Au N ^ Ng +U Nc 



i Ng gcwiu^i 4 / 3 + ^ ion ) + ^Xf 2 ) UN -^ 



^N c ,N a UN c ,N 



7 



Under the condition that N g > 4N C + 1, we have for all (j) € C^(T), 

V(k,l)eK* xTZ* s.t. |fc|,|i|<^iV c , J^I Ng {4>)etei = ^lJ, 
so that, i? MjVc JVb is defined on Vn c by the Fourier matrix 

C,™ 5 ■ — FFT.Af . — - — . FFT,iV 



(l«IV c 



21 



FFT,Ar o 



where, by convention, the last term of the right hand side is equal to zero for k = I. 
We also introduce the variational approximation of (16) 

II™ = inf {^ TFW (^), v Nc \v N f = AAj . (18) 

Any minimizer un c to (18) satisfies the elliptic equation 

^n c u Nc1 (19) 



--^-Au Ne +U Nc 



^C T F\u N f /3 u Nc + V ion u Nc + Vfi° N u $ mh u Nc 



for some A« G 



Theorem 3.1 i<br eac« iV c G N, we denote by un c a minimizer to (18) such that 
(un c ,u) L 2 > and, for each N c G N and N g > 4A^ C + 1, we denote by Ujv c ../v g a 
minimizer to (17) such that (ujv c at J w) i 2 > 0. Then for N c large enough, uw and 
u N c .N g o,re unique, and the following estimates hold true 



\\u Nc - 




c s 7v- (m - s+1/2 - £) 


(20) 


IAjv, - A| 






(21) 


^ rTFW rTFW 

# 




C||ujv c - u||^i 

# 


(22) 






Cs jy3/2+(a-l) + j V -m ) 


(23) 


AAr c ,Af g — Aj\r e | 




CN^ 2 N~ m , 


(24) 


1 rTFW rTFW 1 




C*Ar c 3/2 7V-" 1 , 


(25) 



/or a/Z —to + 3/2 < s < to + 1/2 and for some constants 7 > 0, C > and C s > 
independent of N c and N g . 



Remark 1 More complex orbital-free models have been proposed in the recent years 
[9], which are used to perform multimillion atom DFT calculations. Some of these 
models however are not well posed (the energy functional is not bounded below [1]), 
and the others are not well understood from a mathematical point of view. For these 
reasons, we will not deal with those models in this article. 

Proof [of Theorem 3.1] The estimates (20), (21) and (22) originate from arguments 
already introduced in [2]. For brevity, we only recall the main steps of the proof 
and leave the details to the reader. 
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The difference between (16) and the problem dealt with in [2] is the presence of the 
Coulomb term Dr(M 2 , \v | 2 ), for which the following estimates are available: 



0<D r (p 7 p) 


< 


c\\ P \\i v 

# 


for all p g L|(r), 


(26) 


\Dr(uv, uw)\ 


< 


C|MU»| 


\w\\ L% , forall(z;,i i ;)e(L|(r)) 2 , 


(27) 


\D r (p,vw)\ 


< 


C\\p\\li\ 

# 


\v\\ Ll \\w\\ Ll , for all (p,v,w) G (i|(r)) 3 , 


(28) 


II t /Coulomb || 

II v P 


< 




for all p G L|(r) 


(29) 


t /-Coulomb 1 1 

\ \\ H r 


< 


c\\p\U h 


, for all p G Hjf.(T). 


(30) 



Here and in the sequel, C denotes a non-negative constant which may depend on 
r, V lon and J\f, but not on the discretization parameters. 

LetF(t) = C* TF i 5/3 and f(t) = F'(t) = §CW 2/3 - The function F is in C 1 ([0, +oo))n 
C°°((0, +oo)), is strictly convex on [0, +oo), and for all (ti,t 2 ) 6 M + x E + , 

l/(*§)*2 - /(i?)i 2 - 2/'(i 2 )t 2 (t 2 - < gc TF max(^ /3 , 4 /3 ) |i 2 - t x \ 2 . (31) 

The first and second derivatives of E TFW at the unique positive minimizcr u = y^fi 
of (16) are respectively given by 

(E T ™'(u), v) H i Hl =2(H p0 u,v) 

(E TFW "(u)v,w) H -i ffl = 2{H„ov,w) +4D r (uv,uw) +4 / f(\u\ 2 )\u\ 2 vw, 

* ' * Jr 

where we have denoted by 



:-^LA + f(p) + V im + V p ' 



the TFW Hamiltonian associated with the density p. We recall (see [6] and the 
proof of Lemma 2 in [2]) that (i) u G H™ +1/2 ~ e (T) n C|(r) for each e > 0, (ii) 
u > on K 3 , (hi) A is the ground state eigenvalue of H p o and is non-degenerate. 
Using (26), (27) and the fact that /' > on (0, +oo), we can then show (see the 
proof of Lemma 1 in [2]) that there exist (3 > 0, 7 > and M > such that for all 
v G ffi (T), 

< ((ff„c - X)v, v) H -i H i^ < M\\vf H% (32) 
P\HV ^ ((^ TFW "(") - 2A)!;, v) , Hl < M\\v\\ 2 Hl (33) 

# # ' # # 

and for all n G H^(T) such that IHIl^ = W 1/2 and (w, u) L | > 0, 

7ll« - < (( H P° - - «). (« - u ) V 1 A3 • (34) 

# # ' # 

Remarking that 

£ TFW («aO - £ TFW («) = ((ffpo - A)(« Wc - «), (Ujv c - u))^ Hl 

# # 

+\d t {\u Nc \ 2 -\u\ 2 ,\u N f -\u\ 2 ) 

+ ! F(\u N f)-F(\u\ 2 )-f(\u\ 2 )(\u N f-\u\ 2 ) (35) 
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and using (34), the positivity of the bilinear form Dr, and the convexity of the 
function F, we obtain that 



7-TFW 7-TFW pTFW / \ tttTFWc -\ — ■ | i|2 

i n c -I =E ™{u No )-E (u) > 7||ujv c - u\\ H ±. 



For each N c £ N, um = A/" 1 ^ 2 njv c it/||nAr c w,|| i 2 satisfies un c and IImjvJIi, 2 = Af 1 ^ 2 , 

and the sequence (un c )n c £N converges to u in H^ +1 ^ 2 e (r) for each e > 0. As the 
functional E TFW is continuous on H^(T), we have 

\\m e * < l" 1 (/£f W - I TFW ) < I' 1 (E TFW (u Nc ) - E™(u)) 0. 

# AT C — ► oc 

Hence, (un c )n c gn converges to u in i?i(r), and we also have 

A^ 1 ^ / \^u N f+ J f(\u N f)\u N f + J V ion \u N f + D r (\u N f,\uNf) 
1 



N c —>oo 



X. 



r Jr 

, |Vn| 2 + / /(M 2 )M 2 + / v io >| 2 + iMM 2 ,M 2 ) 

^ ./r 



As /(|ujvJ 2 )ujv c + V ion «iv c + V^° u \° 2 mh u Nc is bounded in L 2 # {T), uniformly in 7Y C , 
we deduce from (19) that the sequence (itjv c )jv c eN is bounded in H%(T), hence in 
L°°(T). Now 



A{u Nc -u) = 2C^ 



Kn c f(\u N f)u Nc - f(\u\ 2 )u + V ion (u Nc - u) + 



t /-Coulomb,,, TT-Coulomb., 
V \u Nc \ 2 U Nc-V\u\2 u 



ion„, I T/Coulomb^ 



+ (i-n We )(/(M> + V Mn u + V[y 



-^nA u n c -u)- {X Nc - X)u 



Observing that the right-hand side goes to zero in L 2 ^ (r) when N c goes to infinity, 

we obtain that (itjv c )./v c eN converges to u in #|(r), and therefore in C°/ /2 {T). 
By a simple bootstrap argument, one can see that the convergence also holds in 
iJ™ +1/2_e (r) for each e > 0. The upper bound in (22) is obtained from (35), 
remarking that 



< J^F(\u Na \ 2 )-F(\u\ 2 )-f(\u\ 2 )(\u Na \ 2 ~\u\ 2 ) 
< yC TF ^max(| Ujv J 4 / 3 , \u\ 4 / 3 )\u Nc - u\ 2 



35 



4/3 



and that 



< ttCtf max ||uatJ.l<*> I \\u Nc - u||? 2 
9 yv c eN 



< D r (\u N f-\u\ 2 ,\u N f-\u\ 2 )<C\\\u N f-\u\ 2 \\l 5 



< AC [ max \\u N Jkoo I ||Mjv c -"||r2 
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The uniqueness of un c for N c large enough can then be checked as follows. First, 
(u n c , Atv c ) satisfies the variational equation 

yv Nc G V Nc , ((H\ u i2 - X N Ju Nc ,v No ) H -i H i = 0. 

# ' # 

Therefore Ajv c is the variational approximation in V/v c of some eigenvalue of H\ UN p. 
As (tiAr c )Ar c6 N converges to it in L°°(r), i?| UJV p — H p o converges to in operator 
norm. Consequently, the n th eigenvalue of H\ UN 12 converges to the n th eigenvalue 
of H p o when N c goes to infinity, the convergence being uniform in n. Together with 
the fact that the sequence (Xn^n^en converges to A, the non-degenerate ground 
state eigenvalue of H p o , this implies that for N c large enough, Ajy c is the ground 
state eigenvalue of i?| UjVo p in Vn c and for all vn c G Vjv c such that \\vn c = A/" 1 / 2 
and (v Nc ,u Nc ) L 2 > 0, 

# 

E TFW (v Nc ) - S TFW (?iAr c ) = ({H\ UN 12 - AjvJ(vjv c - UJV C ), («jv c - u Nc )) H -i Hl 

c # 5 # 

+ ^r(KJ 2 - KJ 2 , |^J 2 - KJ 2 ) 

2\ cVL. |2\ |2\/|„, |2 |„, 1 2 \ 



F(KJ 2 ) - F{\u Ne \') - /(| UW J 2 )(|^J 2 - \u N f) 

> t;\\ v n c - u Nc \\ 2 H i ■ (36) 
z # 

It easily follows that for N c large enough, (18) has a unique minimizcr un c such 

that {um c1 u) L 2 > 0. 

# 

Let us now establish the rates of convergence of |Ajv c — A| and |jujv c — u\\h^- First, 



\n c — A = N 
with 



-1 



{(H\ u \2 - X)(u Nc - u), (u Nc - u)) H -i H i + / w Nc ( u n c - u) 



... „ /(KJ 2 )-/(M 2 ) L , ,2 .^Coulomb/, , , oA 
W N C ~ UN - U ' C ' l u «cl 2 \ UN - +U >- 

We know that the sequence (uAr c )jv e eN converges to u in 7J^ +1 ^ 2 ~ £ (r) and that u > 
in R 3 . Consequently, for N c large enough, the function un c (which is continuous 
and TvL-periodic) is bounded away from 0, uniformly in N c . As / € C°°((0, +00)), the 
function wn c is uniformly bounded in 3 ^ 2 c (r) (at least for N c large enough). 
We therefore obtain that for all < r < m — 3/2, there exists a constant C r G R+ 
such that for all N c large enough, 

|Aat c - A| < C r (\\u Nc - u\\ 2 H1 + \\u Nc - u\\ H -r) . (37) 

In order to evaluate the if^-norm of the error (iijv c — u), we first notice that 

Vwjv G Vzv oJ ||ujv c - < ||«jv c - watJIjji + ||wat c - u\\ H i , (38) 

and that 

\\un c -v Nc \\ 2 h1 < f3- 1 ((E TFW "(u)-2\)(u Nc -v Nc ),(u Nc -v Nc )) H -i Hl 

# # ' # 

= /T 1 ( ((E TFW "(u) ~ 2X)(u Nc - u), (u Nc ~ v Nc )) H -i Hl 

+ ((£ TFW "( W ) - 2A)(u - v Nc ), (u Nc - v Nc )) H -i Hl ) . (39) 

# ' # y 



ii 



For all zn c £ Vn c , 

((E TFW "(u) - 2X)(u Na - u), z Nc ) H -> Hl 

# # 

= -2 Jtf{\u Na \ 2 )u Na - f(\u\ 2 )u - 2f'(\u\ 2 )\u\ 2 (u Nc - u)]z Nc 

-2D r ((u Nc - u)(u Nc + u), (u Nc - u)z Nc ) - 2D r ({u Nc - u) 2 ,uz N J 
+2(X Nc -X) J u Nc z Nc . (40) 

On the other hand, we have for all vm c £ Vn c such that ||wat c ||l 2 = A/" 1 / 2 , 

# 

1, 



u n S u n c - v N J = M - I u Nc v Nc = -\\u No - v N J\ L 2 . 
r Jr A * 

Using (28), (31), (37) with r = and the above equality, we therefore obtain for all 
vn c £ Vn c such that ||«jv ||i 2 =M 1 ^ 2 , 



((E itw (u) ~ 2X)(u Nc - u), (u Na - v Na )) H -, Hl 

# ' # 

< c(\\ u n c -uWh^ \\u Ne -v Ne \\ H ^ 

+ (\\un c -u\\ 2 h i # + \Wn c - u||i|) \\u Nc -«iv c |li|Y (41) 
Therefore, for N c large enough, we have for all vn c £ Vn c such that HvjvJIl 2 = 

\\u Nc - v Nc \\ H i # < C (\\u Nc - uf H i # + \\v Nc - u\\ H iJj . 

Together with (38), this shows that there exists N £ N and C £ R+ such that for 
all N c > N , 

Vv Nc £ Vn c s.t. \\v Nc \\ L 2 = W 1/2 , ||uat c -u\\ H i < C\\v Na -u\\ H i. 

# # # 

By a classical argument (see e.g. the proof of Theorem 1 in [2]), we deduce from 
(2) and the above inequality that 

\\ UNc - U \U<C mm \\v Nc -u\\ H r KCN-l™' 1 ' 2 -*, (42) 
* «N c eVjv c # 

for some constants C independent of N c . 

For w £ i^(r), we denote by tpw the unique solution to the adjoint problem 
[ find tp w £ u such that 



Vv £ u x , ((E TFW "(u) - 2X)iP w ,v) H -i Hl = (w,v) H -i Hl , 



(43) 



where 



u 1 - = l.v£Hl(T) I J tt« = o|. 
The function tp w is solution to the elliptic equation 

-^AV^ + {V ion + vS oulomh + f{u 2 ) + 2f'(u 2 )u 2 — X)tp w + 2V.t 

= 2 1 / f'(u 2 )u 3 ip w + D T (u 2 ,u%l) w ) I u + w - (w,u) L 2 u, 



it 



r 



# 
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from which we deduce that if w £ H^(T) for some < r < m — 3/2, then ip w £ 
H r + 2 {Y) and 

||V>JU+= < CVHIffJ, (44) 
for some constant C r independent of w. Let be the orthogonal projection, for 
the L 2 ^ inner product, of un c on the affine space |w € £#(r) | J r uv =J\f\. One 
has 

from which we infer that 

\ U N C ~ u\\ 2 L 2 = / {u Nc - u)(u* Nc -u)+ {u Nc - u)(u Nc - u* N J 
Jr Jr 



# 



1 



{u Nc - u)(u* Nc - u) + -^jj\\u Na - u\\ L 2 (A/ - - / u Nc u 



2 

r z - /v \ 

(uat c - u)(u* Nc - u) + -y^\\uN c - "llli 

( UjVc - u,u* Nc - u) H -i iH i # + ^\\u Nc - u||£. 



((■E TFW "(u) - 2A)i/; UNc _„,u^ -u) H -i^ + -^\\ u n c - u||l,2 

1 



(( £TFW " - 2A)(uw c ~ m)^™^-^}^- 1 ,^! + TlrH^ ^ "ilk 



{(E TFW "(u) -2\)(u Nc -u),tp UNc - u ) H -i^ +TT7ll«JV c ~ U \\lI 



-jf\\ u Nc-n L i 



V 



AN 

f'(u 2 )u 3 ip UNa ^ u + D r {u 2 ,uip UNa _ u ) 



For all ipN c G Vm c , it therefore holds 

\\ u n c -u\\ 2 L 2 = ({E TFW "(u) - 2X)(u Na - u),tp u u - i'N a ) H -i H i 

# ' # 

+ <(£ TFW "( U ) - 2A)( UNc - u),^ ) H - Sffi + ^\\u Nc u\\ 4 Ll 



+Jf\\ u N c 



f(u 2 )u 3 i/) UNo - u + D r (u 2 ,mp UNc - u ) 

r 



(45) 



Using (28), (31), (37) with r = and (40), we obtain that for all ip Nc £ V Nc H u 1 - , 
((£ TFW ( u )-2A)( WWc - U ),7/w c V lHl < a( ||t*jv — «|||ri 

# ' # y # 

+ ll«iv e - u\\l^ (\\un + hNa - «||z|) WjvJffJ- (46) 



Let us denote by 11^ nu± the orthogonal projector on Vn c H u for the Hi inner 
product and by ip% a = ILy N nuX ipu Nc -u- Noticing that 

Wn.WhI < Uu No -u\\h^ <p- l M\\u No -u\\ Ll , 
13 



we obtain from (33), (45) and (46) that there exists TV € N and C E R+ such that 
for all TV C > TV, 

\\ U N C ~ < C [\\ u n c _ - u\\ L 2 j \\u Nc - + \\u Nc - u\\ H x\\il?v Ne -u ~ tPn c \\h^J ■ 

Lastly, for all v E u 1 - and all TV C G N* 

II* -iW^ < (i+ ^£ Jr J ll-n^||^, (47) 

so that, in view of (2) and (44) 

\\i'u Na -u - iPn c \\h^ < C \\^u Nc -u -njV c ^„ Nc -„||Hi 

< CN- l U UN ^ u \\ Hl 

< CN- x \\u Ne - U \\ Ll . 

Therefore, 

By means of the inverse inequality 

2^ {r ~ s) 



Vvn c £V Nc , \\v Nc \\ H r # < I —j N c r - s \\v Nc \\ H s # , (48) 

which holds true for all s < r and all TV C > 1, we obtain that 

\\u Nc - u\\ H ° # < C s N-( m - s+1 / 2 - e) for all < s < m + 1/2. (49) 

To complete the first part of the proof of Theorem 3.1, we still have to compute the 
ii_^ r -norm of the error {un c — u) for < r < m — 3/2. Let w E H 7 ^(T). Proceeding 
as above we obtain 

w(u Nc -u) = ((E TFW "(u) ~2\){u Nc -i^n^n^^tf-i^i 

+ ((E TFW "(u) - 2X)(u Nc - u),^ w - K Ha nu^ W ) Hi \Hl 

+^||«JV.-u|li| J^f'(u 2 )u 3 ip w +Dr(u 2 ,ui; w ) 

_ W" UW<= - U W 2 Llj r UW - ( 50 ) 

Combining (33), (44), (46), (47), (49) and (50), we obtain that there exists a con- 
stant C E R+ such that for all TV C large enough and all w E H^(T), 

J w{u Nc -u) < C (||ujv e - +N c -( r+1) \\u Nc -u\\ H ^ \\w\\ Hi 
< C N- {m+r+1 ^\\w\\ Hi . 

Therefore 

/ w(u Nc - u) 

\\un c ' u\\ H -r = sup ^— - < CTV-(™+'- +1 /2-), (51) 

* wdHUT)\{Q} \\w\\H r # 
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for some constant C E M+ independent of N c . Using (37), (42) and (51), we end 
up with 

|Ajv-A| < CN-P" 1 - 1 -*). 



Let us now turn to the pseudospectral approximation (17) of (16). First, we notice 
that 



< EY g W (M 1/2 \T\- 1/2 ) 

< c TF AA 5 / 3 |rr 2 / 3 + ||y ion || L o=AA, 

from which we infer that ujv.a^ is uniformly bounded in H^(T). We then see that 



2 



\Vu Na , Ng \'+ / l N AV lon \u NctNg \' + f(\u No , Ng \ 2 )\u Na , Ng \') 
r Jr 



+D r (\u NaiNg \ 2 , |uAr ciA r g | 2 ) 

Using (6), (11) and (26), we obtain that An,n<. also is uniformly bounded. Now, 

Au Nc>Ng = 2C^U No (l Ng (f(\u NetN f)u No>Ng )) + 2C^U No (l Ng (V ion u Nc , Ng )) 
+2C^IL Nc (vfi£ffiu Ne , Ng ) - 2C^\ NciNg u NctNg , (52) 

and we deduce from (4), (6) and (8) that 

\\U Nc (l Ng {f(\u Nc , Ng \ 2 )u Nc , Ng ))\\ L% < {l Ng (f(\u Nc , Ng \ 2 ))) 2 | UN c ,N g \ I 

E (^) 3 f(\ U N c ,N g (x)\ 2 ) 2 \ U N c .N g (x)f 



xGO Ng nr 

5„ „ ,,i/3 / ^ / L x 

3 "- I ^' lLoa I ^ V AT, 



5 



i.zeSjv„nr 



^C T F\\u Nc ,N g || L oo \\uN c .N g || L 4 , 



and that 



\\IL Nc {l Ng (V™u Nc , Ng )) \\ L% < \\Il 2Nc (l Ng (V ioa u Nc , Ng )) \\ L% 

1/2 

< ( / i Ng (\vn 2 \u Nc ,N g \ 2 ) 



r 



< ||U ion ||L~AA 1/2 . 

Besides, using (29), 

ii n *c (vffi%$u Ne , Nt ) w L% < ii^^^j^ 



< Af 1/2 \\u No , Ng \\ 2 
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As UN c ,N g is uniformly bounded in H^(T), and therefore in L%.(T), we get 



\\ u n c mJhI 



\\ u N c ,N g \\ 2 L 2 + \\Au Nc , Ng \\ 2 L 2 



1/2 



< c (i + \\u Ne , Ng \\yi) 

< c(l + \\u NctNg \\^) 



Therefore UN ci N g is uniformly bounded in H^(T), hence in L c 

Returning to (52) and using (9) and a bootstrap argument, we conclude that w_/v c ,jv g 
is in fact uniformly bounded in iJ^ 2+e (T). 

Next, using (36), 

7j\\u Nc ,N g - UN a \\%i 







E7 g W (u Nc ,N g )-E7 g W (u N J 




+ l^(l-2 Ng )(V))(\u Nc , N f 


-WNf) 


+ ^(l-X Ns )(F(| UNciWs | 2 )- 


F(\u N f)) 


J^(l-l Ng )(V))(\u Nc , N f- 




+ |^l-l Ng )(F(\u Nc , Ng \ 2 )- 


F(\u N f)) 



< 



Let g(t, t') = r(l >,~_ t y ' . For N c large enough, uyv c is uniformly bounded away from 
zero; besides, both un c and u^ c pf g are uniformly bounded in H^ 2+e (T) . Therefore, 
g{uN c ,UN c ,N g ) is uniformly bounded in H^ 2+e (T). This implies that the Fourier 
coefficients of ,g(u7v c , ujv c ,jv 3 ) go to zero faster that |fc| -7 / 2 , which implies, using (5) 
and (10), that 

(l-l Ng )(F(\u Nc , N f)-F(\u N f)) 

(1-TjvJ (g(uN c ,U Nc ,N g )) {u Nc ,N a ~ U Nc ) 

< \\U Nc ((l-l Ng )(g(u Nc ,u NciNg )))\\ L2 || u Na , Ng - u Nc \\ L 2 

< CN?/ 2 N-y 2 \\u Nc , Ng - u Nc \\ L% . (53) 
On the other hand, 

((l-l Ng )(V))(\u N ^ Ng \ 2 -\u N f) 

< ||n 2 JV c ((l - ZN g )(V))\\ L 2j\u NciNg + U Nc \\ L ^\\u NciNg -UN c \\lI 

<CN^ 2 N- m \\u Nc , Ng -u Nc \\ Ll . 

Therefore, 

\\u Nc ,N g -u Nc \\ H ^ < CN 3 J 2 N- 7 / 2 . (54) 
We then deduce from (54) and the inverse inequality (48) that (uN c ,N g )N c ,N g >4N c +i 



converges to u in Hi(T), and therefore in L° 



It follows that for N c large 



1G 



enough, UN c .N g is bounded away from zero, which, together with (52). implies that 

(uN a ,N g )N e ,N g >±N e +i is boundei 
therefore be improved, yielding 



(uN c .N g )N c ,N g >4N c +i is bounded in H^ +1 ^ 2 £ (T). The estimates (53) and (54) can 



(l-lN g )(F(\u No , Ng \ 2 )-F(\u N f)) 



r 



< CN 3 J 2 N-( m+1 ^\\u Nc , Ns - u Ne h 



and 

\\u Ne ,N 3 -u Nc \\ H i # <CNV 2 N- m . 

We deduce (23) from the inverse inequality (48). For N c large enough, u Na ^ Ng is 
bounded away from zero, so that f{\uN c ^ g | 2 ) is uniformly bounded in H^ +1 ^ 2 ~ e (T). 
Therefore, the fc th Fourier coefficient of (V ion + f(\u Nc , Ng | 2 )) is bounded by C\k\~ m 
where the constant C docs not depend on N c and N g . Using the equality 



-f r (l-lN g )(V i0 » + f(\u Nc , N /))\ UNc , Ng \ 2 

+Dr(\u Nc , Ng \ 2 , \u Nc , N f - \u N f) + [ (f(\u Nc , N f) ~ f(\u N f))\u Nc , N f 



r 



(23) and (28), we obtain (24). A similar calculation leads to (25). 
Lastly, we have for all vn c G Vn , 

E™(v Nc )-E™»(u Nc , Ng ) (55) 

= (( H u Na . Ng - >^N c ,N g )(vN c ~ U Nc! N ), (v Nc - U Nc . Ng )) H -i H i 

+ 2 D r(|«JvJ 2 - \u Nc ,N g \ 2 , \vnJ 2 - \u Nc ,N g \ 2 ) 

+ E (jf) ( F ^ VN ^ - F d^)| 2 ) - /(I^ (X)| 2 )(|^ (X)| 2 - | UJ V C (X)| 2 )) 

> ((H UNe , Ng ~ ^N c ,N g )(v Nc - u Nc , N ), (v Nc - u NctN )) H -i H i. (56) 

# # 

As Un C! n 9 converges to u in H^(T), the operator N | 2 ~ converges to zero 

in operator norm. Reasoning as in the proof of the uniqueness of un c , we obtain 
that for N c large enough and N g > 4N C + 1, we have for all Vn c € Vn c such that 
IbivJi 2 = A/" 1/2 and {v Nc ,u Nc ) L 2 > 0, 

{(Hu NciNg ~ >^N c ,N g ){vNc ~ UN c ,N g ), (VN C ~ u N c ,N g )) H -\ H ^ > ^|K'JV C - "Nc.jvj^i. 

Thus the uniqueness of un c ,n for iV c large enough. □ 
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